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Abstract 

Effects of an external magnetic field on various properties of the quantum chromodynamics (QCD) 
matter under extreme conditions of temperature and density (chemical potential) have been analysed. 
To this end, we use SU(3) Polyakov linear cj-model and assume that the external magnetic field {eB) 
adds some restrictions to the quarks energy due to the existence of free charges in the plasma phase. 
In doing this, we apply the Landau theory of quantization, which assumes that the cyclotron orbits 
of charged particles in magnetic field should be quantized. This requires an additional temperature 
to drive the system through the chiral phase-transition. Accordingly, the dependence of the critical 
temperature of chiral and confinement phase-transitions on the magnetic field is characterized. Based 
on this, we have studied the thermal evolution of thermodynamic quantities (energy density and trace 
anomaly) and the first four higher-order moment of particle multiplicity. Having all these calculations, 
we have studied the effects of the magnetic field on the chiral phase-transition. We found that both 
critical temperature Tc and critical chemical potential increase with increasing the magnetic field, eB. 
Last but not least, the magnetic effects of the thermal evolution of four scalar and four pseudoscalar 
meson states are studied. We concluded that the meson masses decrease as the temperature increases 
till Tc- Then, the vacuum effect becomes dominant and rapidly increases with the temperature T. At 
low T, the scalar meson masses normalized to the lowest Matsubara frequency rapidly decrease as T 
increases. Then, starting from Tc, we find that the thermal dependence almost vanishes. Furthermore, 
the meson masses increase with increasing magnetic field. This gives characteristic phase diagram of 
T vs. external magnetic field eB. At high T, we find that the masses of almost all meson states 
become temperature independent. It is worthwhile to highlight that the various meson states likely 
have different critical temperatures. 

PACS numbers: 12.39.Fe, 12.38.Aw, 52.55.-s 

Keywords: Chiral Lagrangian, Quark confinement. Magnetic confinement and equilibrium 
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I. INTRODUCTION 

It is believed that at high temperatures and densities there should be phase transition(s) 
between confined nuclear matter and the quark-gluon plasma (QGP), where quarks and gluons 
are no longer confined inside hadron bags |l|. Various theoretical studies have been devoted 
to tackle the possible change in properties of the strongly interacting matter, when the phase 
transition(s) between hadronic and partonic phases takes place under the effect of an external 

5 ? 
















magnetic field {2 7|. It is conjectured that the strongly interacting system (hadronic or partonic) 
can response to the external magnetic held with magnetization, M, and magnetic susceptibility, 
Xm l8|. Both quantities characterize the magnetic properties of the system of interest. Thus, the 
effects of the external magnetic held on the chiral condensates should be rehected in the chiral 
phase-transition jo]. Also, the ehects on the deconhnement order-parameter (Polyakov-loop) 
which includes the conhnement-deconhnement phase-transition can be studied 9|. 

In an external magnetic held, the hadronic and partonic states are investigated in diherent 
models, such as the hadron resonance gas (HRG) model 10| , and other ehective models [U 


The Nambu - Jeno-Lasinio (NJL) model 2( 
quark model 2^ and certain limits of QCD 


22| . the chiral perturbation theory 


23 


19|. 


the 


27| are also implemented. Furthermore, there are 
some studies devoted to the magnetic ehects on the dynamical quark masse s 12811 . The chiral 
magnetic-ehect was studied in context of the Polyakov NJL (PNJL) model j29|. Recently, it 
was reported about lattice QCD calculations in an external magnetic held j^. 


30 


33|. The 


Polyakov linear a-mode’ 
on the system 0,06. 


(PLSM) was implemented to estimate the ehects of the magnetic held 


In the present work, we add some restrictions to the quarks energy due to the existence of free 
charges in the plasma phase. To this end, we apply the Landau theory (Landau quantization) 


36l |. which quantizes of the cyclotron orbits of the charged particles in the magnetic helds. 


We notice that this proposed conhguration requires an additional temperature to drive the 
system through the chiral phase-transition. Accordingly, we hnd that the value of the chiral 
condensates increase with increasing the external magnetic held 0. A few remarks are now 
in order. In many different calculations for the thermal behavior of the chiral condensates 
and the deconhnement order-parameter (Polyakov-loop) using PNJL or NJL j2-4|, the external 
magnetic held was not constant. Also, the dependence of the critical temperatures of chiral 
and conhnement phase-transitions on the magnetic held was analysed Almost the same 
study was conducted in PLSM js-?]. All these studies lead to almost the same pattern, the 
critical temperature of chiral phase-transition increases with increasing the external magnetic 
held. But, the critical temperature of the conhnement phase-transition behaves, OTpositely. 
The latter behavior agrees - to some extend - with the lattice QCD calculations js]. In the 
present work, we study the ehects of external magnetic held on the phase transition and deduce 
the phase-diagram curve using SU(3) PLSM 


In light of this, we recall that the PLSM is widely implemented in diherent frameworks 


and diherent purposes. The LSM was introduced by Cell-Mann and Levy in 1960 39| long 
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time before QCD was known to be the theory of strong interaction. Many studies have been 


performed with LSM like 0{A) LSM 


39j, 0(4) LSM at hni 


U{Nf)i LSM for Nj = 2, 3 or even 4 quark flavors 42l-l45[| . In order to obtain reliable results, 


;e temperature 


40|,l4lj andU{Nf)rX 


Polyakov-loop corrections have been added to LSM, in which information about the conhning 
glue sector of the theory was included in form of Polyakov-loop potential. This potential is 


to be extracted from the pure Yang-Mills lattice simulations 46 


491]. So far, many studies 


were devoted to investigating the p 


Polyakov-loop forms with two M, 


rase diagram and the thermodynamics of PLSM at different 


5l| and three quark flavors 


38 


52 


53| . Also, the magnetic 


held effect on t 
PLSM 


3, 


le QCD phase-transition and other system properties are investigated using 


3J,l35|. 


The present paper is organized as follows. In section HIl we introduce details about SU(3) 
PLSM under the effects of an external magnetic held. Section IIIII gives some features of the 
PLSM in an external magnetic held, such as the quark condensates, Polyakov loop, some 
thermal quantities, the phase-transition(s) and scalar and pseudoscalar meson masses under 
the magnetic held ehect. In section HVl the hnal conclusions and outlook shall be presented. 


II. APPROACH 


The Lagrangian of LSM with Nf = 2 + 1 quark havors and Nc = 3 color degrees of freedom. 


where the quarks couple to the Polyakov-loop dynamics, was introduced in Ref. 38|, [52, l53[ | . 


C = Cchiral - W(0, 0*,T), 


( 1 ) 


where the chiral part of the Lagrangian Cchirai = has SU{3)l x SU{3)ji symmetry 


54j . l55| . The Lagrangian with Nf = 2 -|- 1 consists of two parts. The hrst part represents 
fermions, Eq. (Ej) with a havor-blind Yukawa coupling g of the quarks. The coupling between 
the ehective gluon held and quarks, and between the magnetic held, B, and the quarks is 
implemented through the covariant derivative 


^9 = - gTa{aa + i'y57Ta))i>f, 


( 2 ) 


where the summation runs over the three havors, / = 1, 2, 3 for u-, d- and s-quark, respec¬ 
tively, Ta is the Gell-Man matrices. The havor-blind Yukawa coupling, g^ should couple the 
quarks to the mesons |56j. The coupling of the quarks to the Euclidean gauge held, A^, was 
















































discussed in Ref l47| . For the Abelian gauge field, the influence of the external magnetic 
held, is given by the covariant derivative [^, 


-Dm = d^-iA^-iQ A 


EM 
M ’ 


(3) 


where = 5 ^ A“A“/2 and A^^ = (0, Bx, 0, 0) and Q is a matrix dehned by the quark electric 
charges Q = diag(g^i, qa, qs) for up, down and strange quarks, respectively. The interaction of 
charged pion tt^ = (tti ± ^ 2 ) /\/2 with the magnetic held is included by = d^j, — ie A^ with 
e is the electric charge 

The second part of chiral Lagrangian stands for the the mesonic contribution, Eq. l]T|), 


Cm = Tr(a^<hta^$ - - Ai[Tr(<Ft$)]2 

- A 2 Tr($t $)2 + c[Det($) + Det(<f)^)] + Tr[Tr($ + <h^)]. 


(4) 


In Eq. (jl]), <F is a complex 3x3 matrix, which depends on the a a and 


55| . where 7 ^ are 


the chiral spinors, Ua are the scalar mesons and vTa are the pseudoscalar mesons. 

The second term in Eq. ([I]), W(0, 0*, T), represents the Polyakov-loop ehective potential [46(], 


which is expressed by using the dynamics of the thermal expectation value of a color traced 
Wilson loop in the temporal direction <F(T) = {V{x))/Nc. Then, the Polyakov-loop potential 
and its conjugate read (j) = (TtcV)/Nc, and (j)* = (TrcD^j/W,, respectively. V, which stands 
for the Polyakov loop, can be represented by a matrix in the color space 


V{x) = Vexp 


i / dTA4{x,T) 


(5) 


where 13 = 1/T is the inverse temperature and A 4 = iA^ is the Polyakov gauge 
Polyakov loop matrix can be given as a diagonal representation 571. 


46 


47|. The 


In the PLSM Lagrangian, Eq. CDi, the coupling between the Polyakov loop and the quarks 
is given by the covariant derivative of — iA^ It is apparent that the PLSM 

Lagrangian is invariant under the chiral flavor-group. This is similar to the original QCD 
Lagrangian SSHGOj. In order to reproduce the thermodynamic behavior of the Polyakov loop 
for pure gauge, we use a temperature-dependent potential f/(0, 0*,T). This should agree with 



present work, we use U{(j),(j)*,T), Landau-Ginzburg type potential, as a polynomial expansion 



























(6) 


in (j) and (jf SsMgI] 


W(0 ,r,T) h{T) 




T4 2 ^ 

where, <p and cf)* are introduced previously and b 2 (T) = ao + ai (To/T) + a 2 (To/T)^ + a 3 (Tq/T)^, 
where constants are qq = 6.75, oi = —1.95, 02 = 2.625, 03 = —7.44, 63 = 0.75 and 64 = 7.5. 

In Eq. dnD, the Vandermonde Jacobian contribution, k ln[J(0, 0*)], was ignored due the 
small value of k. In principle, the Vandermonde term comes from the change of variables from 
vector potential to 0 in the path integral and should guarantee a reasonable behavior of the 


mean field approximation 


72j, i.e. it was suggested to solve the problem that the normalized 


Polyakov loop becomes greater than 1 at very high temperatures. 

^ [1 - 600 *+ 4(03 + 0 * 3 ) 

where J(0, 0*) is the Vandermonde determinant, which is not explicitly space-time dependent. 
The dimensionless parameter k would be dependent on the temperature and the chemical 
potential. Therefore, k should be estimated, phenomenologically. 

In order to reproduce the pure gauge QCD thermodynamics and the behavior of the Polyakov 
pop as a function of temperature, we use the parameters listed out above in this section 


59l |. In calculating the grand potential, we use the mean field approximation 
fi(T,/i) = V(cr„,crj^)+W(0,0*,T) + fi^^(T;0,0*,5). 
The purely mesonic potential is given as. 


(7) 


m 


U (da,, ay) = —{ai + ay) - h^a^, - hya. 




2^2 


+ + o (2-^1 + + t(Ai + A2)cry, 


( 8 ) 


62| and based on Landau 


where m^, hy, Ai, A 2 and c are the model fixed parameters [5^. The quarks and anti 
quark contribution to the medium potential was introduced in Ref 
quantization and magnetic catalysis concepts, App. we get 

\qf\BT^ fdp 


^ u=0 '' 


In 

In 


1 + 3(0 + 0*e 
1 + 3(0* + 
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It is worthwhile to highlight that the chemical potential used everywhere in the manuscript is 
the quark one, fij with / being the quark flavor. The different variables are elaborated in the 
App. 1^ The potential at vanishing eB reads 

{In [1 + 3{(j) + 

+ In [1 + 3(0* + } . (10) 

This is the system free of Landau quantization. 

The Landau theory quantizes of the cyclotron orbits of charged particles in magnetic held. 
For small magnetic fields, the number of occupied Landau levels (LL) is large and the quantiza¬ 
tion effects are washed out, while for large magnetic fields, the Landau levels are less occupied 
and the chiral symmetry restoration occurs for smaller values of the chemical potential. 
According to Eqs. (jO]) and (fTOj) . Eq. ((^ get an additional term. 


n{T,Hf,eB) = U{a^,ay) +U{(j),(j)*,T) + n^^{T,Hf;(j),(j)*,eB) + (11) 


where fif] 0, 0*) represents the potential term at vanishing magnetic held, (5o,es switches 

between the two systems; one at vanishing and one at finite magnetic held. 

We notice that the sum in Eqs. ([6]), ([9]) and ([8]) give the thermodynamic potential density 
as in Eq. (JTj). By using the minimization condition, App. |Bl we can evaluate the parameters. 
Having the thermodynamic potential, Eq. ([71), we can determine all thermal quantities including 
the higher-order moments of particle multiplicity, and then mapping out the chiral phase- 


diagram 


The meson masses are defined by the second derivative with respect to the 


corresponding fields of the grand potential, Eq. ([7|), evaluated at its minimum. 


III. RESULTS 


The results of the chiral condensates and ay, section [III A1 the thermodynamic quantities, 
section IIIIB1 the non-normalized and normalized higher-order moment of particle multiplicity, 
section IIII Cl and section IIII C 21 respectively, the chiral phase-transition, section IIII PI and 
Anally the meson masses, section IIII El are introduced as follows. 
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A. Phase transition: qnark condensates and order parameters 


The thermal evolution of the chiral condensates, and ay, and the Polyakov order param¬ 
eters, (j) and (j)* is calculated from Eq. ([7]) at finite chemical potential and finite magnetic held 
using the minimization conditions given in Eq. flBip . The dependence on the four parame¬ 
ters, temperature T, chemical potential p, magnetic held B and minimization parameter with 
respect to it the minimization condition shall be analysed. 

In left-hand panel (a) of Fig. [H the normalized chiral condensates, ax and ay, are given 
as function of temperature at vanishing chemical potential and diherent magnetic held values, 
eB = 10 MeV^ (double-dotted curve), 200 MeV^ (solid curve) and 400 MeV^ (dotted curve). We 
notice that both condensates increase with increasing the magnetic held, eB. This dependence 
seems to explain the increase in the chiral critical temperature with the magnetic held. This 


- in turn - agrees with various studies using PLSM and PNJL 
moderated (smoother) with increasing magnetic held. 


SB 


The condensates become 
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Fig. 1: (Color online) Left-hand panel (a): the normalized chiral-condensates, ax (lower curves) and 
ay (upper curves), are given as function of temperature at vanishing chemical potential and different 
magnetic field values, eB = 10 MeV^ (double-dotted curve), 200 MeV^ (solid curve) and 400 MeV^ 
(dotted curve). Right-hand panel (b): the same as in left-hand panel but at a constant magnetic field 
eB = 200 MeV^ and different quark chemical potentials, p = 100 MeV (solid curve), 200 MeV (dotted 
curve) and 300 MeV (double-dotted curve). 




The right-hand panel (b) of Fig. [T] shows the chiral condensates, ax and ay, as function 
of temperature at constant magnetic held eB = 200 MeV^, and hnite chemical potentials, 
p = 100 MeV (solid curve), 200 MeV (dotted curve) and 300 MeV (double-dotted curve). 
Both condensates decrease with increasing the chemical potentials. This dependence gives a 
signature for the decreasing behavior of the chiral critical temperature with increasing the 


Q 










chemical potential, which obviously agrees with our previous calculations 
become rowdy (sharper) with increasing chemical potential. 


. The condensates 



Fig. 2: (Color online) Left-hand panel (a): the Polyakov-loop field and its conjugation, (p and p*, are 
given as function of temperature at a vanishing constant chemical potential and different magnetic 
field values, eB = 10 MeV^ (double-dotted curve), 0.2 GeV^ (solid curve) and 400 MeV^ (dotted 
curve). Right-hand panel (b): the same as in the left-hand panel but at a constant magnetic field, 
eB = 0.2 GeV^ and different quark chemical potential values, p = 0.1 GeV (solid curve), 0.2 GeV 
(dashed curve) and 0.3 GeV (double-dotted curve). 


The left-hand panel (a) of Fig. [2] shows the Polyakov-loop held and it is conjugation, 
(p (upper curves) and (p* (lower curves), as function of temperature at a vanishing chemical 
potential and different magnetic held values, eB = 0.01 GeV^ (double-dotted curve), 0.2 GeV^ 
(solid curve) and 0.4 GeV^ (dotted curve). Both helds decrease with increasing the magnetic 
held. This behavior explains the dependence of the conhnement critical temperature on the 
magnetic held. At vanishing chemical potential, (p = (p*. Both Polyakov-loop helds become 
smoother with increasing magnetic held. 

The right-hand panel (b) draws the same as in left-hand panel but at a constant magnetic 
held eB = 0.2 GeV^ and diherent quark chemical potential values, /i = 0.1 GeV (solid curve), 
0.2 GeV (dashed curve) and 0.3 GeV (double-dotted curve). We hnd that (p increases with 
increasing the chemical potential values but <p* decreases. This behavior seems to agree with 
our previous calculations 38(]. At hnite chemical potential, (p > (p*. 


We conclude that the Polyakov-loop helds, (p and 0*, increase with T, Fig. [2l At vanishing 


/i, both <p and (p* decrease with increasing eB. At hnite /i, we hnd that <p increases, while cp'' 
decreases with eB. 






B. Thermodynamic quantities 


In this section, we introduce some thermal quantities like energy density and trace anomaly. 


As we discussed in Ref. 


, the purely mesonic potential, Eq. 


gets inhnity at very low 


temperature and entirely vanishes at high temperature. From this numerical estimation, we 
concluded that this part of potential is only effective at very low temperatures. Its dependence 
on the external magnetic held has been checked and was found that hnite eB comes up with 
very tiny contribution to this potential part. As the present study is performed at temperatures 
around the critical one, this potential part can be removed from the effective potentials given 
in Eq. ([7]). In Eq. ([8]), the chiral condensates, a’s, are small at hnite temperature. Fig. 
[U Therefore, much smaller values are expected for their higher orders and multiplications. 
Opposite situation is likely at very small temperatures. 

1. Energy density 


The energy density, e/T'^, at hnite quark chemical potential, /x/, can be obtained as 


e(r,/i/,eR) = 


d 


In Z(T, jjf, eB). 


( 12 ) 


S(i/r) 

In section IIII A1 we have estimated the parameters, the two chiral condensates, ax and ay and 
the two order parameters of the Polyakov-loop and it’s conjugation, cj) and 0*, respectively. 
Thus, we can substitute all these into Eq. (fT2|) . 

The left-hand panel (a) of Fig. [3]presents the normalized energy density, e/T^, as function of 
temperature at vanishing chemical potential. In calculating the results, Eqs. ([9]) and Eq. (na 
are implemented as given in Eq. flTTD . The general tenmerature-dependence is not absent. Also, 
we notice that e/T^ is sensitive to the change in eB 62|]. Increasing eB seems to increase the 
critical temperature, at which the system undergoes phase transition. As the chiral condensates 
become smoother with increasing eB, the thermodynamic quantities, such as energy density, 
behave accordingly, i.e. the phase transition becomes smoother as well. 

The right-hand panel (b) shows e/T^ as function of temperature at a constant magnetic 
held eB = 0.2 GeV^ and varying quark chemical potentials, /r = 0.1 GeV (long-dashed curve), 
0.2 GeV (dash-dotted curve) and 0.3 GeV (double-dotted curve). The solid curve represents 
the results in absence of an external mag netic held but at jj, = 0.1 GeV. We note that e/T^ 


is not as sensitive to the change in /r 38| as to the external magnetic held. Despite the lack 
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Fig. 3; (Color online) Left-hand panel (a): the normalized energy density, ejT"^, is given as function 
of temperature at a vanishing chemical potential and different magnetic fields, without magnetic field 


(solid curve) |38|], eB = 0.01 GeV^ (long-dashed curve), eB = 0.2 GeV^ (dashed curve) and eB = 
0.4 GeV^ (double-dotted curve). The right-hand panel (b) shows the same as in the left-hand panel 
but at a constant magnetic field eB = 0.2 GeV^ and different chemical potential values, ^ = 0.1 GeV 
(long-dashed curve), 0.2 GeV (dash-dotted curve) and 0.3 GeV (double-dotted curve). The upper 
curves represent results from Eq. [9] plus Eq. [TOl while lower curves are based on thermodynamic 
derivatives from the thermal potential, Eq. O 


of chemical potential dependency, which can be understood due to the large magnetic field 
applied, it is believed to affect contrary to the chemical potential. To this indirect dependency 
of /i and eB, we shall devote a separate work. Again, it seems that increasing Hf, decreasing 
T 

C* 


2. Trace anomaly 


At hnite quark chemical potential, the trace anomaly known as interaction measure reads 

(13) 


- 3p(T,iif,eB) ^ p(T,iif,eB) 

2^4 ffip 2^4 


In Fig. m we notice that the normalized trace-anomaly under the effect of an external magnetic 


held becomes smaller than the corresponding quantity in absence of magnetic held 


at high 


temperature. This can be explained due the restrictions added to the quark energy by the 
Landau quantization through the magnetic held. We hnd that increasing eB increases the 
critical temperature. This behavior can be understood because of the dependence of the chiral 
condensates. Fig. [T]and the Polyakov-loop potential. Fig. [2] on eB. 

In the left-hand panel (a), the trace anomaly (e —3p)/T^, is given as function of temperature 


19 





















T at vanishing chemical potential but different values of the magnetic fields, vanishing 
(solid curve), eB = 0.01 GeV^ (long-dashed curve), eB = 0.2 GeV^ (dash-dotted curve) and 
eB = 0.4 GeV^ (double-dotted curve). We notice that the trace anomaly increases with T until 
the chiral symmetry is restored. Then, increasing T reduces the normalized trace anomaly. 
The peak represents the critical temperature corresponding to a certain magnetic field. We 
find that Tc increases with increasing eB. 


w 




(a) eB=0.00 GeVj 

eB=0.01GeV, 
H - 0.00 GeV eB=0.20 GeVj 
eB=0.40 GeV^ 








0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.^ 
T[GeV] 



Fig. 4: (Color online) The same as in Fig. [3] but for the the trace anomaly (e — 3p)/T^. 


The right-hand panel (b) of Fig. 0] shows the same as in the left-hand panel but at a constant 
magnetic field eB = 0.2 GeV^ and different chemical potentials, /i = 0.0 GeV (solid curve), 
0.1 GeV (long-dashed curve), 200 MeV (dash-dotted curve) and 300 GeV (double-dotted curve). 
The trace anomaly increases with T until the chiral symmetry is fully restored. The peaks are 
positioned at of the certain value for chemical potential. Here, we find that decreases 
with increasing /i. The sensitivity to /i is not as strong as to eB. This might be interpreted 
as the high magnetic field applied seems to contradict the effects of the chemical potential. 
In other words, should the magnetic field adds energy to the system, the chemical potential 
requires energy in order to produce new particles. We notice that the dependence on the quark 
chemical potential is more obvious that that shown in Fig. |3l 


3. Magnetic catalysis effect 


In App. lAl we discuss the magnetic catalysis, Eq. flAbll . and the so-called dimension 
reduction concepts, Eq. flAbll . Due to the effects of the magnetic field, the latter would mean 
modifying the sum over the three-dimensional momentum space to a one-dimensional one. 
According to Ref. |38|, the effect of this reduction reduces also the value of the quantity by 


almost two third from the expected value. This would explain the difference between results 

















at vanishing and that at hnite eB, left-hand panels (a) of Figs. [3] and IH for instance. In 
the present work, we distinguish between two types of systems. In the first one, the Landau 
quantization should be implemented, i.e. taking into account the magnetic effects, while in 
the other system, the external magnetic held is not taken into consideration, i.e. no magnetic 
contribution to the thermal system. 


C. Higher-order moment of particle multiplicity 


38 


641 as 


The higher-order moment of the particle multiplicity is dehned 

a /i* T4 ’ ^ ^ 

where the pressure p(T,fi,B) = —TO lnZ{T, p, B)/dV is related to the partition function, 
which in tern is related to the potential. In Z{T, fi, B) = —V Q{T, fi, B)/T. 

In this section, we introduce the hrst four non-normalized moments of the particle multiplic¬ 
ity calculated in PLSM under the effects of an external magnetic held. The thermal evolution 
is studied at a constant chemical potential but diherent magnetic helds and also at a constant 
magnetic held but diherent chemical potentials. Doing this, it is possible to map out the chiral 
phase-diagram, for which we determine the irregular behavior in the higher-order moments as 
function of T and p. 


1. Non-normalized higher-order moments 


Here, we introduce the non-normalized higher-order moments of the particle multiplicity 


38l |. The left-hand panels (a) of Figs. [5], [6], [3 and |8] show the hrst four non-normalized 


moments of the quark distributions. These quantities are given as function of temperature at a 
constant chemical potential p = 0.1 GeV and diherent magnetic helds, eB = 0.1 GeV^ (double- 
dotted curve), 0.4 GeV^ (dashed curve) and 0.7 GeV^ (dotted curve). We hnd that increasing 
temperature rapidly increases the four moments. Furthermore, the thermal dependence is 
obviously enhanced, when moving from lower to higher orders. The values of the moment are 
increasing as we increase the magnetic held. The huctuation in the third- and fourth-order 
moments rehect the increase of the critical temperature with increasing the magnetic held. 
The critical temperature can, for instance, be dehned where the peaks are positioned. 

The right-hand panels (b) of Figs. O [6], [7] and [H present the same as in the left-hand 
panels but at a constant magnetic held eB = 0.1 GeV^ and diherent chemical potentials. 
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Fig. 5: (Color online) Left-hand panel (a): non-normalized quark number density, is given as 
function of temperature at a constant chemical potential = 0.1 GeV and different values of the 
magnetic field, eB = 0.1 GeV^ (double-dotted curve), 0.4 GeV^ (dashed curve) and 0.7 GeV^ (dotted 
curve). Right-hand panel (b) shows the same as in the left-hand panel but at a constant magnetic 
field eB = 0.1 GeV^ and different chemical potentials, /i = 0.01 GeV (double-dotted curve), 0.1 GeV 
(dashed curve) and 0.2 GeV (dotted curve). 
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Fig. 6: (Color online) The same as in Fig. [5] but for quark number susceptibility m 2 . 


/i = 0.01 GeV (double-dotted curve), 0.1 GeV (dashed curve) and 0.2 GeV (dotted curve). It 
is apparent that increasing temperature rapidly increases the four moments of quark number 
density. Furthermore, the thermal dependence is obviously enhanced, when moving from lower 
to higher orders. The values of the moment are increasing as we increase the chemical potential. 
But the critical temperature Tc decrease with /r. The peaks are positioned at the critical 
temperature. 
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Fig. 7: (Color online) The same as in Fig. Obut for the third-order moment of quark number density 
m3. 
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Fig. 8: (Color online) The same as in Fig. [5]but for the fourth-order moment of quark number density 

m4. 


2. Normalized higher-order moments 


The statistical normalization of the higher-order moments requires a scaling of the non- 
normalized quantities, section [Til C 1[ with respect to the standard deviation cr, which is related 


to the susceptibility x oi the fluctuations 


63 


64l | in the particle multiplicity. It is conjectured 


that the dynamical phenomena could be indicated by large fluctuations in these dimensionless 


moments and therefore, the chiral phase-transition can be mapped out 


63|. Due to the sophis¬ 


ticated derivations, we restrict the discussion here to dimensionless higher-order moments 


This can be done when the normalization is done with respect to the temperature or chemical 
potential. 

The higher-order moments of the particle multiplicity normalized with respect to tempera¬ 
ture are studied in dependence on the temperature at a constant chemical potential and different 
magnetic fields. Also they are studied at different chemical potentials and a constant magnetic 
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field. The corresponding expressions were deduced in Ref. 
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Fig. 9: (Color online) Left-hand panel (a): the dimensionless quark number density is given 

as function of temperature at constant chemical potential fx = 0.1 GeV and different magnetic fields, 
eB = 0.1 GeV^ (double-dotted curve), eB = 0.4 GeV^ (dashed curve) and eB = 0.7 GeV^ (dotted 
curve). Right-hand panel (b) shows the same as in the left-hand panel but at a constant magnetic 
field eB = 0.1 GeV^ and different chemical potential values, // = 0.01 GeV (double-dotted curve), 
fi = 0.1 GeV (dashed curve) and fx = 0.2 GeV (dotted curve). 
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Fig. 10: (Color online) The same as in Fig. [9] but for the dimensionless quark number susceptibility 

In left-hand panel (a) of Figs. [HI [lU] and [TT] the hrst three normalized moments are given 
as function of temperature at a constant chemical potential /i = 0.1 GeV and different mag¬ 
netic fields, eB = 0.1 GeV^ (double-dotted curve), 0.4 GeV^ (dashed curve) and 0.7 GeV^ 
(dotted curve). We hnd that the values of the moments are increasing as the magnetic held 
increases. The huctuations in the normalized moments would dehne the dependence of the 
critical temperature on the magnetic held. 

The right-hand panels (b) of Figs. [9], [10] and [11] present the hrst three normalized moments 
as function of temperature but at a constant magnetic held eB = 0.1 GeV^ and diherent 
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Fig. 11: (Color online) The same as in Fig. Obut for the dimensionless third-order moment of the 
quark number density m^/T. 


chemical potentials, /r = 0.01 GeV (double-dotted curve), 0.1 GeV (dashed curve) and 0.2 GeV 
(dotted curve). We notice that the moments of quark multiplicity increase with the chemical 
potentials. That the peaks at corresponding critical temperatures can be used to map out the 
chiral phase-diagrams, T vs. eB and T vs. /i. 


D. Chiral phase-transition 


Now we can study the effects of the magnetic field on the chiral phase-transition. In a 


previous work 


, we have introduced and summarized different methods to calculate the 


critical temperature and chemical potential, Hc, by using the fluctuations in the normalized 
higher-order moments of the quark multiplicity or by using the order parameters. The latter is 
implemented in the present work. The PLSM has two order-parameters. The first one presents 
the chiral phase-transition. This is related to strange and non-strange chiral condensates, 
and ay. The second one gives hints for the confinement-deconfinement phase-transition, the 
Polyakov-loop fields, 0 and <p*. Therefore, for the models having Polyakov-loop potential, we 
can follow a procedure as follows. We start with a constant value of the magnetic field. By using 
strange and non-strange chiral-condensates, a dimensionless quantity reflecting the difference 
between the non-strange and strange condensates Ag ,,(T) as a function of temperature at 
fixed chemical potentials will be implemented. This procedure give one point in the T-/i- 
chart, at which the chiral phase transition takes place. At the same chemical potential as in 
previous step, we deduce the other order-parameter related to the Polyakov-loop fields as a 
function of temperature. These calculations give another point {in T and /i chart), at which 
the deconfinement phase-transition takes place. By varying the chemical potential, we repeat 


1 8 















these steps. Then, we hnd a region (or point), in (at) which the two order-parameters, chiral 
and deconhnement, cross each other, i.e. equal each other. It is assumed that such a point 
represents phase transition(s) at the given chemical potential. In doing this, we get a set of 
points in a two-dimensional chart, the QCD phase-diagram. 



Fig. 12; (Color online) The chiral phase-diagram, T/Tco vs. /i/pcO) at different values of the magnetic 
fields, eB = 0.1, 0.2, 0.3, 0.4 and 0.5 GeV^ from top to bottom. The normalization quantities 
Tco = 0.15 MeV and pco = 0.3 GeV were deduced in Ref. 381. 


In Fig. [121 we compare five chiral phase-diagrams, T/T^o vs. /i//ico, with each others at 
eB = 0.1, 0.2, 0.3, 0.4 and 0.5 GeV^ from top to bottom. T/Tco are plotted against p/pco. 


where t 
Ref. 


le two normalization quantities Tco = 0.15 GeV and pco = 0.3 GeV were deduced from 
This should give an indication about the behavior of the critical temperature and the 


critical chemical potential of the system under the effect of the magnetic field. Apparently, we 
conclude that both critical temperature and critical chemical potential increase with increasing 
the magnetic field. 


E. Meson masses 

The masses can be deduced from the second derivative of the grand potential, Eq. ([7]) with 
respect to the corresponding fields, evaluated at its minimum, which is estimated at vanishing 










expectation values of all scalar and pseudoscalar fields 
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(15) 


where a and h range from 0, • • • ,8 and and are scalar and pseudoscalar mesonic fields, 
respectively. Obviously, i stands for scalar and pseudoscalar mesons. 

The scalar meson masses 651 
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with {ml)ab = m? 5ah - QGabc(^c + 9.J=abcd(yc(yd- The expressions for Qabc and Tabcd can be 


found in Ref. 


65| . On the tree level, can be determined according to dU{a)/daa = 0 = 


nr ^Qabc^b T ^ 2 bcd ^d 9^a' 
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Fig. 13: (Color online) Left-hand panel (a): the scalar meson masses, m„ from Eq. (|16p . nif^ from Eq. 
(I17p are given as function of temperature at a constant magnetic field eB = 0.1 GeV^ and different 
chemical potentials, p = 0.0 GeV (dotted curve), 0.1 GeV (dashed curve) and 0.2 GeV (double-dotted 
curve). Right-hand panel (b): the same as in left-hand panel (a) but for from Eq. (ITSP and 

niag from Eq. (unp. 


In Fig. [131 fhe scalar meson masses, from Eq. (IT6|) . and nif^ from Eq. (ITT)) are given 
as function of temperature at a constant magnetic field eB = 0.1 GeV^ and different chemical 














potentials /i = 0.0 GeV (dotted curve), 0.1 GeV (dashed curve) and 0.2 GeV (double-dotted 
curve). We conclude that the scalar meson masses decrease as the temperature increases. This 
remains until T reaches the critical value. Then, the vacuum effect becomes dominant and 
rapidly increases with the temperature. The effect of the chemical potential is very obvious. 
The masses decrease with the increase in chemical potential. This explains the phase diagram of 
temperatures and chemical potentials at a certain magnetic field. The decrease of the critical 
temperature with increasing chemical potential is represented by the bottoms (minima) in 
thermal behavior of meson masses before switching on the vacuum effect. 

In Fig. [m the four scalar meson masses, rUa from Eq. flT^ . m/p from Eq. flTT)) . m^-^g 
from Eq. fll8p and m^g from Eq. flT^ are given as function of temperature at two values of 
chemical potential, /i = 0.1 GeV left-hand panel (a) and p = 0.2 GeV right-hand panel (b) and 
different magnetic fields, eB = 0.1, 0.2, 0.3, 0.4, 0.5 and 0.6 GeV^ from top to bottom. We 
notice that the scalar meson masses decrease as the temperature increases, until it reaches the 
critical temperature. Then, the vacuum effect gets dominant and apparently increases with the 
temperature. The effect of magnetic field is very obvious. The masses increase as the magnetic 
field increases. This explains the phase diagram of temperatures and magnetic field at a certain 
chemical potential. 

In Fig. [15], the normalized scalar meson masses, mg- from Eq. (IT6H . m/p from Eq. (IT7|) . 
from Eq. flT8|) and mo-g from Eq. flT^ are given as function of temperature at two values 
of chemical potential, p = 0.1 GeV in left-hand panel (a) and /i = 0.2 GeV in right-hand 
panel (b) and different magnetic fields, eB = 0.1, 0.2, 0.3, 0.4, 0.5 and 0.6 GeV^ from top to 
bottom. The normalization is done due to the lowest Matsubara frequencies, 2ttT, App. O At 
high temperatures, we notice that the masses of almost all meson states become temperature 
independent, i.e. constructing a kind of a universal bundle. This would be seen as a signature 
for meson dissociation into quarks. In other words, the meson states undergo deconfimement 
phase-transition. It is worthwhile to highlight that the various meson states likely have different 
critical temperatures. 

At low temperatures, the scalar meson masses normalized to the lowest Matsubara frequency 
rapidly decrease as the temperature increases. Then, starting from the critical temperature, we 
find that the thermal dependence almost vanishes. The magnetic field effect is clear, namely 
the meson masses increase with increasing magnetic field. This characterizes T vs. eB phase 
diagram. 
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Fig. 14: (Color online) Left-hand panel (a): the four scalar meson masses are given as function of 
temperature at a constant chemical potential /r = 0.1 GeV and different magnetic fields, eB = 0.1, 
0.2, 0.3, 0.4, 0.5 and 0.6 GeV^ from top to bottom. Right-hand panel (b): the same as in left-hand 
panel but at chemical potential /r = 0.2 GeV. 
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Fig. 15: (Color online) Left-hand panel (a): the scalar meson masses normalized with respect to the 
lowest Matsubara frequency are given as function of temperature at a constant chemical potential 
fj, = 0.1 GeV and different magnetic fields, eB = 0.1, 0.2, 0.3, 0.4, 0.5 and 0.6 GeV^ from top to 
bottom. Right-hand panel (b): the same as in left-hand panel but at chemical potential = 0.2 GeV. 






































The pseudoscalar meson masses [65 1 
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Fig. 16: (Color online) The same as in Fig. [13] but for pseudoscalar meson masses, from Eq. (1201) . 

rrip from Eq. m, rurj^g from Eq. (f22l) and rupg from Eq. (|23ll . 


In Fig. [T6l the pseudoscalar meson masses, from Eq. fl20l) . from Eq. fl2Tl) . from 
Eq. fl22l) and m^g from Eq. fl2^ are given as function of temperature at a constant magnetic 
held eR = 0.1 GeV^ and different chemical potentials, p = 0.0 GeV (dotted curve), 0.1 GeV 
(dashed curve) and 0.2 GeV (double-dotted curve). It is obvious that the pseudoscalar meson 
masses remain constant at low temperature. At temperatures > Tc, the vacuum effect becomes 
dominant. Accordingly, the pseudoscalar meson masses increase with the temperature. We 
shall notice that even contribution by the vacuum will be moderated through the normalization 
with respect to the lowest Matsubara frequency. 

In Fig. [T71 the four pseudoscalar meson masses, from Eq. fl20l) . from Eq. fl2Tl) . 
from Eq. (|22|) and from Eq. (|23|) are given as function of temperature at two constant 
chemical potentials, p = 0.1 GeV in left-hand panel (a) and p = 0.2 GeV in right-hand panel 














2.5 


^ 1-5 

I 

1 

0.5 

0 


eB-O.lGeV^ . 

3 

eB-O.lGeV^ . 

H-U.lGeV jB^o.2GeV^ . 


^1 = 0.2 GeV ,B,o.2GeV^ . 

eB=0.3GeV^ . 


eB=0.3GeV^ . 

eB=0.4 GeV“ 


eB=0.4 GeV’ 

eB=0.5GeV‘ 


eB=0.5GeV“ • • 

cR-nf^rw- . 

a 


_ _ 

§ 1.5 


'. ' 

1 

1 



0.5 

n . 

, 

0 

, 


T[GeV] 



T[GeV] 


(b) 

^ = 0.2 GeV 


eB=0.1 GeV^ 
eB=0.2 GeV“ 
eB=0.3 GeV- 
eB=0.4 GeV“ 
eB=0.5 GeV“ 
eB=0.6 GeV- 







eB-O.l GeV^ " 



eB=0.2GeV^ ■■ 



eB=0.3GeV- 



eB=0.4 GeV" 



eB=0.5GeV' 



eB=0.6GeV^ 





ns 







0.2 0.3 

T[GeV] 


‘ (b) 

eB-0.1 GeV^ -• 


p = 0.2 GeV 

eB=0.2GeV- 



eB=0.3GeV- . 



eB=0.4 GeV" 



eB=0.5GeV" 



eB=0.6GeV^ • 





ns 




0.2 0.3 

T[GeV] 



3 


B = 0.1 GeV :b=o:2gV . 

2J 

:b=o:2gV . 

eB=0.3GeV^ . . 


eB=0.3GeV^ . . 

eB=0.4 GeV^ 


eB=0.4 GeV^ 

eB=0.5 GeV^ . 


eB=0.5GeV^ - - 

eB-0.6 GeV^ . 


eB-0.6 GeV^ " 


“ 15 






1 

.“. 

Ins 

1 

nNS . 


0.5 

0 



0.2 0.3 

T[GeV] 


0.2 0.3 

T[GeV] 


Fig. 17: (Color online) Left-hand panel (a): the pseudoscalar meson masses are given as function of 
temperature at a constant chemical potential fi = 0.1 GeV and different magnetic fields, eB = 0.1, 0.2, 
0.3, 0.4, 0.5 and 0.6 GeV^ from top to bottom. Right-hand panel (b) shows the same as in left-hand 
panel but at a constant chemical potential /r = 0.2 GeV. 


(b) and different magnetic fields, eB = 0.1, 0.2, 0.3, 0.4, 0.5 and 0.6 GeV^ from top to bottom. 
Again, at low temperature, the masses remain temperature-independent. At T > Tc, the 
vacuum effect is switched on. Accordingly, the masses increase rapidly with the temperature. 

In Fig. [THl the four pseudoscalar meson masses, from Eq. (|20|) . from Eq. dZD 
from Eq. and from Eq. fl2^ normalized with respect to the lowest Matsubara fre¬ 
quency are given as function of temperature at two constant chemical potentials, /r = 0.1 GeV 
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Fig. 18: (Color online) Left-hand panel (a): the pseudoscalar meson masses normalized to the lowest 
Matsubara frequency are given as function of temperature at a constant chemical potential /r = 
0.1 GeV and different magnetic fields, eB = 0.1, 0.2, 0.3, 0.4, 0.5 and 0.6 GeV^ from top to bottom. 
Right-hand panel (b): shows the same as in left-hand panel but at a constant chemical potential 
^ = 0.2 GeV. 


in left-hand panel (a) and fi = 0.2 GeV in right-hand panel (b) and different magnetic helds, 
eB = 0.1, 0.2, 0.3, 0.4, 0.5 and 0.6 GeV^ from top to bottom. The normalization should result 
in temperature-independent mesonic states. This would be seen as a signature for meson disso¬ 
ciation into quarks. It is obvious that various critical temperatures can be assigned to various 
pseudoscalar meson states. The normalized masses starts with high values reflecting conhne- 
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ment, especially at low temperatures. Then, they decrease as the temperature increases until 
the critical temperature, Tc, which differs for different meson states. At higher temperatures, 
the dependence of meson masses on temperature is almost entirely removed. 


Symbol 

PDG [73] 

Present Work 

PNJL [74, la 

Lattice QGD 

Hot QGD [76] 

PAGS-GS [77] 

TT 

134.970 ±6.9 

120 

126 

134 ±6 

135.4 ±6.2 

K 

497.614 ± 24.8 

509 

490 

422.6 ± 11.3 

498 ± 22 

rj 

547.853 ± 27.4 

553 

505 

579 ± 7.3 

688 ± 32 

/ 

V 

957.78 ± 60 

965 

949 

- 

- 


aetween pseudoscalar meson masses, = 0 ■*■, in MeV and the corresponding 


Tab. I: A comparison_ 
results from PNJL 
Group (PDG) and the lattice QGD calculations 


75l |. Both are compared with the ex 


aerimental measurements, the Particle Data 

r^. 


As introduced in Ref. 78|], Tab. [Ipresents a comparison between pseudoscalar meson nonets 
in various effective models, like PLSM (present work) and PNJL [7^ confronted to the particle 
data group [73[ and lattice QCD calculations [76|, l77| . The comparison for scalar states would 
be only partly possible. Some remarks are now in order. The errors are deduced from the 
htting for the parameters used in calculating the equation of states and other thermodynamics 
quantities. The output results are very precise for some of the lightest hadron resonances 
described by the present model, PLSM. An extended comparison is given in Ref. 781. 


IV. CONCLUSIONS AND OUTLOOK 


The QCD phase-diagram at vanishing chemical potential and hnite temperature subject 
to an external magnetic held gained prominence among high-energy physicists, for instance. 


our previous work 


62| was based on two concepts in order to explain the effects of external 


magnetic held on the QCD phase-diagram. Another study was done in the framework of Nambu 
- Jona-Lasinio (NJL) model and Polyakov NJL (PNJL) model 6^. The main idea is that the 
scalar coupling parameter is taken dependent on the magnetic held intensity. Thus, the scalar 
coupling parameter decreases with the magnetic held increase, we also implemented the relation 
between the magnetic held and the scalar coupling parameter in order to ht for the lattice QCD 


results 


3l| . We conclude the increase in the magnetic held increases the critical temperature. 

































In the presence work, we use the Polyakov linear cr-model and assume that the external 
magnetic field adds some restrictions to the quarks due to the existence of free charges in the 
plasma phase. In doing this, we apply Landau theory (Landau quantization), which quantizes 
of the cyclotron orbits of charged particles in magnetic field. First, we have calculated and 
then analysed the thermal evolution of the chiral condensates and the deconfinement order- 
parameters. We notice that the Landau quantization requires additional temperature to drive 
the system through the chiral phase-transition. Accordingly, we find that the value of the 
chiral condensates increase with increasing the external magnetic field. In the contrary to 
various previous studies, the effects of the external magnetic field are analysed, systematically. 
Accordingly, the dependence of the critical temperatures of chiral and confinement phase- 
transitions on the magnetic field could be characterized. We deduced T-/i curves using SU(3) 
PLSM in external magnetic field. 

Furthermore, by using mean field approximation, we constructed the partition functions and 
then driven various thermodynamic quantities, like energy density and interaction rate (trace 
anomaly). Their dependence on temperature and chemical potential recalls to highlight that the 
effects of external magnetic field on the chemical potential was disregarded in all calculations 
at finite chemical potential. 

We have analysed the first four non-normalized higher-order moments of particle multiplicity. 
The thermal evolution was studied at a constant chemical potential but different magnetic fields 
and also at a constant magnetic field but different chemical potentials. Doing this, the chiral 
phase-diagram can be mapped out. We determined the irregular behavior as function of T and 
/i. We found that increasing temperature rapidly increases the four moments and the thermal 
dependence is obviously enhanced, when moving from lower to higher orders. The values of the 
moment are increasing as we increase the chemical potential. But the critical temperature Tc 
decrease with /i. The peaks are positioned at the critical temperature. 

The higher-order moments normalized to temperature are studied at a constant chemical 
potential and different magnetic fields. Also they are studied at different chemical potentials 
and a constant magnetic field. The statistical normalization requires scaling with respect to the 
standard deviation, a, where a is related to the susceptibility y or the fluctuations. Due to the 
sophisticated derivations, the discussion was limited to dimensionless higher-order moments. 
This can be done when the normalization is done with respect to the temperature or chemical 
potential. We find that the higher-order moments increase with the magnetic field. We found 
that the moments increase with the chemical potentials. That the peaks at corresponding 


critical temperatures can be used to map out the chiral phase-diagrams, T vs. eB and T vs. /i. 

The effects of the magnetic field on the chiral phase-transition have been evaluated. There 
are different methods proposed to calculate the critical temperature and chemical potential, 
fic, through implementing fluctuations in the normalized higher-order moments or by the order 
parameters. The latter was implemented in the present work. It is obvious that PLSM has 
two types of order-parameter. The first one gives the chiral phase-transition and is related 
to strange and non-strange chiral condensates. The second one gives hints for deconfinement 
phase-transition. Therefore, we can follow a procedure that at a constant magnetic field and 
by using strange and non-strange chiral-condensates, a dimensionless quantity would reflect the 
difference between the non-strange and strange condensates Ag^s{T) as a function of tempera¬ 
ture at fixed chemical potentials, i.e. chiral phase-transition. At the same chemical potential, 
we can deduce the other order-parameter related to the Polyakov-loop fields as function in tem¬ 
perature. Both calculations give one point, at which the two order-parameters crossing each 
other. It is assumed that such a point represents the transition point at the given chemical 
potential. We repeat this at various chemical potentials and get a set of points in a two- 
dimensional chart, the QCD phase-diagram. We have compared hve QCD phase-diagrams, 
T/Tco vs. /i//ico, with each others at hve different values of the magnetic held. We found that 
both critical temperature and critical chemical potential increase with increasing the magnetic 
held. 

The masses can be deduced from the second derivative of the grand potential with respect 
to the corresponding helds, evaluated at its minimum, which is estimated at vanishing expec¬ 
tation values of all scalar and pseudoscalar helds. We have studied scalar and pseudoscalar 
meson masses as function of temperature at two diherent values of magnetic held and diher- 
ent chemical potentials. We concluded that the meson masses decrease as the temperature 
increases. This remains until T reaches the critical value. Then, the vacuum ehect becomes 
dominant and rapidly increases with the temperature. The decrease of the critical temperature 
with increasing chemical potential is represented by the bottoms (minima) in thermal behavior 
of meson masses before switching on the vacuum ehect. At low temperatures, the scalar me¬ 
son masses normalized to the lowest Matsubara frequency rapidly decrease as the temperature 
increases. Then, starting from the critical temperature, we hnd that the thermal dependence 
almost vanishes. Furthermore, the meson masses increase with increasing magnetic held. This 
characterizes T vs. eB phase diagram. At high temperatures, we notice that the masses of 
almost all meson states become temperature independent, i.e. constructing kind of a universal 


line. This would be seen as a signature for meson dissociation into quarks. In other words, the 
meson states undergo deconfimement phase-transition. It is worthwhile to highlight that the 
various meson states likely have different critical temperatures. 


Appendix A: Magnetic catalysis 


For simplicity, we assume t 
From the magnetic catalysis 


ha 


the direction of the magnetic field B goes along 2 ;-direction. 


67l | and by using Landau quantization, we find that when the 


system is affected by a strong magnetic field, the quark dispersion relation will be modified to 
be quantized by Landau quantum number, n > 0, and therefore the concept of dimensional 
reduction will be applied. 


Eu = 

yjpl + ml + |g„|(2n + 1 - a)B, 

(Al) 

Ed = 

^Jpl + ml+ \qd\{2n + l - a)B, 

(A2) 

Es = 

y/pl + 771"^+ \qs\{2n -f 1 - a)B, 

(A3) 


where a is related to the spin quantum number and S {a = ±S/2). Here, we replace 2n + l — a 
by one quantum number z/, where i/ = 0 is the Lowest Landau Level (LLL) and the Maximum 


Landau Level (MLL) was determined according to Eq. flA7p 68| . ruf, where / runs over u-, d- 
and s-quark mass. 


rriq = g- 


2 ’ 


rris = g 




(A4) 

(AS) 


We apply another magnetic catalysis property j67|, namely the dimensional reduction. As 
the name says, the dimensions will be reduced as D —> D — 2. In this situation, the three- 
momentum integral will transformed into a one-momentum integral 

^ \qf\BT f dp 

v=0 


T 


2-k 


(2-lM- 


(A6) 


(27r)^ 27r 

when 2 — represents the degenerate in the Landau level, since for LLL we have single 
degenerate and doublet for the upper Landau levels, 

_ ^^QCD , . 

“ 2\qj\B- ^ ^ 

We use niq and for non-strange and strange quark mass, i.e. the masses of light quarks 

degenerate. This is not the case for the electric charges. In section [TTl qu, qd and are 

elaborated. 


















Appendix B: Minimization condition 


We notice that the thermodynamic potential density as given in Eq. 0, which has seven 
parameters m^,hx,hy,Xi,X 2 ,c and g, two unknown condensates ax and ay and the order pa¬ 
rameters for the deconhnement, 0 and 0*. The six parameters hx,hy, Xi, X 2 and c are hxed 


in the vacuum by six experimentally known quantities 


55| . In order to evaluate the unknown 


parameters ax, ay, 0 and 0*, we minimize the thermodynamic potential, Eq. ([7]), with respect 
to ax, ay, 0 and 0* or 0ij and 0|j. Doing this, we obtain a set of four equations of motion, 

dVti dVti (9Di 


= 0 , 


dax day 00 00 * 

meaning that ax = (fx, ay = dy, (j) = (j) and 0* = 0* are the global minimum. 


(Bl) 


Appendix C: Matsubara frequencies 


In hnite temperature held theory, the Matsubara frequencies are a summation over the 
discrete imaginary frequency, Sy = T where g{iujn) is a rational function, ojn = 

2n7rT for bosons and ujn = (2n-|- l)7rT for fermions and n = 0, ±1, ±2, • • • is an integer 
playing the role of a quantum number. By using Matsubara weighting function hy{z), which 
has simple poles exactly located a.t z = ioJn, then 

T 


5'r, = ^ ^ g{z) hy{z) dz. 


(Cl) 


where p = ± stands for the statistic sign for bosons and fermions, respectively. hy{z) can be 
chosen depending on which half plane the convergence is to be controlled. 


hyi^Z') 


= 


l+nr,{z) 

7 X ’ 


'I T ’ 


(C2) 


where ny{z) = (l -|- is the single-particle distribution function. 

The mesonic masses are conjectured to have contributions from Matsubara frequencies 


ra|. 


Furthermore, at high temperatures, > T^, the behavior of the thermodynamic quantities, in¬ 
cluding the quark susceptibilities, the masses is affected by the interplay between the lowest 
Matsubara frequency and the Polyakov loop-correction 70|. We apply normalization for the 
different mesonic sectors with respect to lowest Matsubara frequency 7l| in order to charac¬ 
terize the dissolving temperature of the mesonic bound states. It is found that the different 















mesonic states have different dissolving temperatures. This would mean that the different 
mesonic states have different Tc’s, at which the bound mesons begin to dissolve into quarks. 
Therefore, the normalized masses should not be different at T > Tc. To a large extend, their 
thermal and dense dependence should be removed, so that the remaining effects are defined by 
the free energy |^, i.e. the masses of free besons are defined by m/. 

That the masses of almost all mesonic states become independent on T, i.e. constructing 
kind of a universal line, this would be seen as a signature for meson dissociation into quarks. 
It is a deconfinement phase-transition, where the quarks behave almost freely. In other words, 
the characteristic temperature should not be universal, as well. So far, we conclude that the 
universal characterizing the QCD phase boundary is indeed an approximative average (over 
various bound states). 
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